Semi-classical formula for quantum tunneling in asymmetric double-well potentials 
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Despite quantum tunneling has been studied since the advent of quantum mechanics, the literature 
appears to contain no simple (textbook) formula for tunneling in generic asymmetric double-well 
potentials. In the regime of strong localization, we derive an succinct analytical formula based on the 
WKB semi-classical approach. Two different examples of asymmetric potentials are discussed: when 
the two localized levels are degenerate or not. For the first case, we also discuss a time-dependent 
problem showing quantum Zeno effect. 
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I. INTRODUCTION 

Quantum tunneling is of continuing interest in many 
contemporary areas of physics [l[ . The simplest problem 
can be formulated as a degree of freedom x whose poten- 
tial energy V(x) has a double- well shape. In the classical 
limit and zero temperature one expects that an initial 
state prepared in one well is stable. Quantum tunneling 
allows the possibility to escape from one side to the other 
passing under the classically forbidden region. Canonical 
examples are an nitrogen atom in ammonia molecule [H 
or an electron in a double quantum dot Q. 

However, in recent years, there has be a break- 
through in the experimental study of macroscopic quan- 
tum tunneling, which is being used to create and 
study " Schrodingcr-cat" states. Examples include Bose- 
Einstein condensates in a double trap [3-01 and quan- 
tum superconducting circuits based on Josephson junc- 
tions |8Hl4|. Macroscopic quantum tunneling is impor- 
tant to test the validity of the quantum mechanics on 
scales larger than the atomic one [l5j]. The investiga- 
tion of these fundamental issues will be also useful for 
advanced technological applications, such as the devel- 
opment of devices for quantum information processing 



These experiments involve a system tunneling from one 
macroscopic state to another. Despite the complexity of 
this process, it is often remarkably well described by the 
physics of a single particle in a double-well potential in 
which the variable x corresponds to a collective macro- 
scopic variable. 

The double- well needs not be symmetric, and in many 
experiments the asymmetry can be changed by modifying 
externally some tunable parameters (see, for instance, 
Refe.jlia- 14|. Ref.Q and references therein). 

One can formulate the problem as a potential Vr,(x) 
whose shape depends on a dimensionless parameter rj 
that quantifies the asymmetry, i.e. rj = corresponding 
to the symmetric case. In the limit of high energy barrier, 
i.e. Vq S> fko where Vb is the energy scale of the barrier 
and w the typical harmonic frequency in the wells, the 
low-energy physics reduces to the standard two-level sys- 
tem: £l(??), sr(t)) are the energies of two localized states 



coupled by quantum tunneling with amplitude v{rj). 

The last parameter 1^(77) has to be determined from 
the given double-well potential. Remarkably, the litera- 
ture appears to contain no simple (textbook) formula for 
tunneling in generic asymmetric double- well potentials 
fnl - liHl ]. Few exceptions are the works in Refs.fl4. [23 - [25j 
in which general methods based on sophisticated tech- 
niques are discussed but useful analytical formulas are 
presented only for specific shapes of the potential. 

In this work we revisit the problem. By using the stan- 
dard WKB-approach, we demonstrate that it is possible 
to express the amplitude v{rj) as the simple formula: 



v (77) = A (77) yJv L {ri)v R {ri) 



A{ri) 



fV(Q)-e L (r))Y , fV(Q)-e R (ri) 



\V{0)-e r (v)J \V(0)-e L (r,) 



(1) 

(2) 



where and v R {ri) are the tunneling amplitudes 

associated to two symmetric double-well potentials: 
Vl{x,T)) = Vl(—x,T}) and Vr(x,t]) = Vr(—x,T)). As 
shown in Fig. [TJ they are defined by the equations 
V v (x) = Vl(x,t]) for x < and V v (x) = Vr(x,t)) for 
x > 0, Fig. [1] V^(0) is the maximum of the potential at 
x = 0. 




Figure 1. (Color on line) The asymmetric potential V(x, rj), 
the solid (black) line, is conceived as the result of merging 
two symmetric double- well potentials: V v (x) = VL{x,r]) for 
x < 0, the dashed (red) line, whereas V v (x) = Vn(x,ri) for 
x > 0, the dotted (blue) line. 

Now and VR{rj) can be easily obtained by using 

the well-known analytic formula for symmetric double- 
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wells (see Refs.[26M28|) reported in Eqs. ([T3J , (JUJ) of this 
article for completeness. 

Hence the analytical and succinct formula Eqs. (JIJ, ([2]) 
together with the Eq. (|13|) allow a direct calculation of 
the tunneling amplitude in an asymmetric double-well. 
To the best of our knowledge, this general formulation 
has never been proposed although it can be of great use 
to experimentalists in quantifying their results. 

The rest of the paper is organized as follows. In 
Sec. HI] wc introduce the two-level system and we recall 
the equivalence of the two semi-classical methods, i.e. 
the WKB approximation and the instanton technique. 
Then, by using the simpler WKB method, wc derive the 
formula Eqs. ([U,© in Sec. Mil As an example of ap- 
plication, wc discuss two cases in Sec. IIVI In the first 
one wc consider an asymmetric double-well potential in 
which the degeneracy is removed as the asymmetry is in- 
troduced: the bias quartic potential with sl{vi) 7^ £r(v) 
for ij ^ (see Fig. [3^). In the second case, we consider 
a particular situation in which the asymmetry of the po- 
tential is introduced without removing the degeneracy 
£h{il) = £r{vi) = £ (sec Fig. [5b). Now the role of the 
asymmetry in the tunneling dynamics appears in a clear- 
cut way as the Rabi frequency is directly related to the 
tunneling amplitude ft£l = 2v(rf). For this second case, 
we also discuss a time-dependent problem corresponding 
to a special case of the quantum Zeno effect {23 . |30| in 
an asymmetric double-well potential in which the asym- 
metry varies with time rj(t). 

The appendix contains comparisons between the exact 
numerical results and the semi-classical formula for the 
case of the two potentials discussed in this paper. 



a) 1 

>° 




-2 -1 1 2 -2-101234 

X / a X / a 



Figure 2. (Color on line) Two examples of asymmetric double- 
well potential at different values of 77: r\ = (symmetric) full 
lines, 77 — 0.5 dashed (red) lines, (a) The quartic potential 
Eq. (J23J). (b) The parabolic potential of Eqs. (p^ . (f27)) . p9|) 
with el — sr — —Vo/2. The horizontal lines represent 
schematically the two energies £1,(77) and en(rj). 



II. MODEL AND APPROXIMATIONS 
A. Shape of the potential 

For the potential V(x), we assume that its maximum is 
set at the origin V(0) = 0, the minima at left x = < 
and at right x = a R > 0, and that the potential increases 
in the limit \x\ — > 00 so that the eigenstates {ip n (x)} 



decay far away the origin and the energy spectrum {E n } 
is discrete. 

Here we focus on the low-energy regime of tunneling. 
Then, a priori, only a few eigenstates are involved in the 
dynamics. Specifically, one can restrict to the two lowest 
energy states Eq,Ei: the ground state tpo (x) and the first 
excited state ipi(x). 

B. Two-level effective model 

Around the two minima s = L,R, whose local har- 
monic frequencies are uj s , it is possible to solve locally 
the Schrodingcr equation to obtain two localized states 
4>l{x) and 4>r(x). They are well localized into the wells 

over a range of order a s = [fi,/(27nw s )] 1 ^ 2 (m is the par- 
ticle's mass) when the heights associated to the left and 
right energy barrier are large compared to the kinetic 
energies of localization 

[1/(0) - V(a s )} > ^ ~ fkj s . (3) 

However 4>l(%) and 4>r{ x ) also decay inside the barrier 
as we explain below. 

Their energies e s = V(a s ) + fuv s /2 arc assumed to be 
close to the energy range spanned by Eq,E\. Wc assume 
that the asymmetry is not too strong to yield resonance 
between the first (approximate) excited state in one well 
and the (approximate) ground state of the opposite well 
[23| . Increasing the asymmetry, the two- level description 
is no more valid as the resonant condition is approached, 
for instance when Sl(v) — £ r(v) + ^flX 7 / - )- 

In the regime defined by Eq. the linear combina- 
tion of the two wave functions 4>l(x) and 4>r(x) is to be 
the best choice as a (real) two-terms expansion to ap- 
proximate the exact ground state and the first excited 
state. The most general linear combination reads 

V>o(ar) = cos(6)cj) L {x) - sin(6)<j> R (x) , (4) 
ifa.(x) = sm(0)<j> L (x) + cos(8)(t) R (x) . (5) 

Eqs. (HI),© provide the condition of orthonormality be- 
tween ipQ and ipi assuming that 4>l{%) and 4>r(x) are 
themselves orthogonal and normalized. The above dis- 
cussion motivates the introduction of the standard two- 
level Hamiltonian 

H= J2 e s (v)\s){s\-v(v)J2\ s )( s '\ > ( 6 ) 

s—L 7 R s^s' 

which describes quantum tunneling between the two lo- 
calized states in the left well and in the right one and v{rf) 
is the amplitude for tunneling. The full Hilbert space is 
thus spanned by only two states. Solving the equation for 
the two eigenstates ipo and ipi of the 2x2 matrix Eq. (|6|) , 
one obtains the relation between the angle 9 appearing 
in Eqs. dU),© and the parameters e s , v 

sin (26) = , / ^ . . (7) 

V ^(r?) + [A £ (r?)/2] 2 
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where Ae(r;) = £1,(77) — £r(t))- For Ae = the rotat- 
ing angle is 9 = 7r/4. Then the general Hamiltonian 
for two-level system with degeneracy Eq. ([6]) has always 
eigenstates corresponding to symmetric or antisymmet- 
ric linear combinations of the two states L and R. It is 
worthwhile to stress that this result holds even if the left 
and right wave functions 4>l(x) and </>_r(x) are generally 
different for an arbitrary asymmetric potential. On the 
other hand, the asymmetry affects the tunneling ampli- 
tude v(r]) as we explain below. 

Given the parameters entering the effective Hamilto- 
nian Eq. (j6|), any physical quantities can be evaluated. 
For instance, solving the time-dependent Schrddingcr 
equation and assuming \ip(t = 0)) = \L) at the initial 
time, we obtain Pn(t), the probability to have the system 
at right and at the time t. With the boundary condition 
-Pr(O) = 0, its derivative satisfies the following equation 



dP R (t) 
dt 



sin(26>) sin (Sit) 



where the Rabi oscillation frequency is hSl = E\ 
and E Q , E x = (Ae( V )/2) T [(Ae( V )/2) 2 + v 2 ( V )f 2 . 



(8) 



have the relation 

dtfi(x) 



d.r 



x=0 



= -fc(O)0(O). 



(11) 



When the expression Eq. (| 10[) is inserted in Eq. we 
obtain the result 



v = —k(0)(j) 2 (0) 



(12) 



As discussed in Ref . [28j . expanding Eqs. (|10[) . ([T2"T) around 
the singular point x, it is possible to recover exactly the 
instanton solution [26|, [23 in which we have the expo- 
nential integrals extending from the first x = — a to the 
second minimum x = +a. For the sake of completeness, 
we recall here the result: 



v = hui\ 



e C e - k I-: dx ^2m[V(x)-V{a)} 



where the numerical prefactor C is 



C 



dx 



V y / 2m[V(x) - 7(a)} a-x 



(13) 



(14) 



C. WKB approximation and instanton technique 
for symmetric double-well 

To obtain the amplitude v for quantum tunneling in 
the semi-classical regime defined by Eq. flHJ, the well- 
known approaches are the instanton technique and the 
WKB approximation. 

The WKB method is based on the Lifshitz-Herring for- 
mula [il], H3 • For a symmetric double- well, it reads 




in which 4>(x) represents the (approximated) localized 
solution in the double- well potential (4> — 4>l — <Pr)- 

In the past the discrepancy of the results given by 
the two methods has been extensively discussed. For 
different potentials, they resulted in a difference by the 
well-known factor (e/n) 1 ^ 2 [33|, |34|- For a symmetric po- 
tential, Garg showed that this difference is not a failure 
of the formula Eq. © but it was related to the choice 
of the wave function <j> approximating the local solution 
[28j. Introducing x, as the classical turning point where 
the local energy e = El = £r intersects the potential 
barrier e = V(x), the corrected WKB solution inside the 
barrier (0 < x < x) for the symmetric double-well reads 

mm 

(t>{x) ^2^) e ~ iIxSdx ' Kx ' ) ' k[x) = ^(H^-e), 

(10) 

where k(x) is now the inverse of the local penetration 
length. As dV(x)/dx = at the maximum at x = 0, we 



D. Discussion 

Despite their equivalence for the formula of the quan- 
tum tunneling amplitude, the instanton technique re- 
mains particularly advantageous when the particle is cou- 
pled to an external environment [36| . Then the quantum 
dissipative dynamics can be formulated in terms of the 
path integral in which the relevant object is the Euclidean 
action of the full system (and not the wave functions) 
[37| . In the presence of coupling with a dissipative exter- 
nal bath, a particle moving in an asymmetric double- well 
potential can even relax towards its minimal configura- 
tion USUI. 

The coupling with the environment is weak as long as 
7a 2 /h <C 1 where 7 is the linear friction coefficient [381 ] . 
When the coherent quantum dynamics of the particle in 
the potential is of interest (the decoherence time is much 
greater than the Rabi frequency), the effective two- level 
model holds and the tunneling amplitude is an intrinsic 
quantity which characterizes the quantum isolated sys- 
tem. 

As Bcnderskii ct al. demonstrated that the two semi- 
classical methods are equivalent even for asymmetric 
double- well potentials [24|], one can choose to use the 
simpler method to tackle the problem, i.e. the WKB 
approach. 

III. SEMI-CLASSICAL FORMULA FOR 
ASYMMETRIC POTENTIALS 

We now formulate the WKB-approach for an arbitrary 
asymmetric potential. The starting point is the conser- 
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vation law for the probability density p(x,t) = \^(x,t)\ 
dp(x,t) dJ(x,t) 



dt dx ' 

in which the probability current reads 

J(x,t) = Im V(x,t) V-^- 

m \ dx 



(15) 



(16) 



Integrating the Eq. (fl"5"j) from x — to x — +oo, at the 
left hand-side we have the derivative of the probability 
to have the particle in the right space x > 0, -Pr(£) 

dt m \ dx J x =o 

(17) 

Assuming ^(x, 0) = 4>l{x) at the initial time, by invert- 
ing Eqs. ((H), (JS|), we have a simple expression for the evo- 
lution of the wave function 

$(x,t) = cos(9)e-i Eot iP o (x)+sm(0)e--k Elt M x ) ■ (18) 
We now insert Eq. (fl~8)) into Eq. (fl7)) and we obtain 
dP R (t) 



dt 



— sm(20) sm(fit) i>o -j— 

2m \ dx dx J X _ Q 



A sin (20) s in(m)U L ^ 



^) (19) 

dx J x=0 



where we have used again Eqs. (HI),© in the last line. 
Comparing Eq. (fH))) with Eq. ^ we obtain the relation 
between the tunneling amplitude and the two wave func- 
tions: 



2m \ dx 



dx 



(20) 



x=0 



For a symmetric potential, the left and right states are 
equal 4>r{x) = 4>{x) = 4>l{—x) and we recover the stan- 
dard formula for the tunneling amplitude Eq. (|9|). Wc 
now introduce x s , the crossing point of the left and 
right energies with the potential barrier V(xl) = £l and 
V(xr) = er. Recalling Eq. (fTU)) . the left and right wave 
functions (s = L, R) inside the barrier are given by 



4> s {x) 



muj s J e -if* Ldx ' kL( - x "> (x<x L ) 
2nek s {x) X 1 e"* f: R d *' k *^ ( x > 



(21) 

1 /2 

where fc s (x) = [2m(7(i) — e 8 )] . Owing to our choice 
of the maximum's position, we have dV(x)/dx = at 
x = so that the derivative dcj> s /dx is related to the 
function itself </> s (0), e.g. Eq. (fTTT) . Then the tunneling 
amplitude Eq. ([20]) reads 



" =tt- (*«(O)0l(O)^h(O) + fcz,(0)^(0)^(0)) 
2 m 




v^(o)Mo)fe(o)^(o) 



\JvrvZ 



(22) 



In the last equation, we have used the formula Eq. (fP2|) 
for the tunneling amplitude i/£ and z/r for two symmet- 
ric potentials defined by Vl(x) = V(x) for x < and 
Vk(a;) = V(x) for a; > 0. The Eq. is the result 
Eqs. (HJ),© as k 8 (0) = y/2m(V(0) - e s ). 



IV. APPLICATIONS 
A. Asymmetric quartic potential 

The first example of asymmetric potential is the quar- 
tic potential including a linear term 



V(x) = V 




9 9 

muj a 



(23) 

The positions of the minima for 77 = are at ±a, the 
harmonic frequency at the bottom of the wells is u. For 
the symmetric case, the tunnel amplitude reads 



v=A[~) y/(Hu)V e 



16 yq_ 

3 hu> 



(24) 



When the asymmetry parameter is introduced 77 > 0, 
a liner term is added in the potential which also re- 
moves the energy degeneracy (see Fig. [2^,). For < 77 < 
8/(3-\/3) the potential is has still a maximum and two 
minima which are now shifted from x = and x = ±a. 
They are given by the formula 



2n 



— arctan - 



,3^3*7 



(25) 

with n = 2 for the maximum ac and n = 1 and n = 
for the left ol and right minimum clr. The harmonic 
frequency are given by (uj s /u>) = (3/2)(a s /a) — 1/2 
with s — L, R. 

After changing the origin of the x axis, x' = x—ac, and 
knowing of the position of the two minima, we can use 
the formulas Eqs. (P),©- The prefactor Eq. © can be 
easily calculated by observing that it corresponds to A = 

(l/2)[(a/C«) 1/4 + (Wa) 1/4 ] with C« = V(a c )-V(a s )- 
hw s /2 for s = L,R. We use the formula for the tunneling 
amplitude in the symmetric double- well Eqs. (pTH|) . (|14[) to 
obtain Vl for and vr from the two potential Vl{x) = 
V(x) for x < a c and Vr(x) = V(x) for x > a c 

The result is given in Fig. [3] where the ratio between 
the tunneling amplitude ^(77) of the asymmetric potential 
and the amplitude for the symmetric case ^(0) is shown 
for different values of Vo/hoj. 

Here we summarize the results. First, we observe that 
the asymmetry of the potential always reduces the tun- 
neling amplitude for the quartic potential. This behav- 
ior is not universal but it depends on the specific shape 
of the potential and how it is deformed (see next sec- 
tion). Second, the resonant condition el = £r + fiujR is 
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Figure 3. (Color on line) The tunneling amplitude v(r\) for 
the asymmetric quartic potential Eq. ()23|) as a function of the 
asymmetry parameter rj scaled with the symmetric value ^(0) 
at different values of the ratio Vo/ftui. 



matched at the values 77 = 0.12,0.25,0.5 for the ratios, 
respectively, Vq/(Hlj) = 4, 2, 1 shown in Fig. [3] This con- 
dition reads Vq/(Hlj) < 1/(277) to the leading order in 77 
and for Vq/(Hw) ^> 1 (see below and the Appendix). Ap- 
proaching these values, the two-level description becomes 
insufficient. Third, the renormalization of the tunneling 
amplitude due to the asymmetry increases when increas- 
ing the ratio Vq/Hlj. This behavior can be explained by 
looking at vl and vr. These two amplitudes have ex- 
pressions similar to the symmetric case v, Eq. ([24| . The 
leading dependence of 77 is in the exponential term, i.e. 
the term S of Eq. (fT5|) which has a scale factor Vq/Huj. 

We conclude by observing the energy difference Ae 
overwhelms the corrections of the asymmetry in the 
tunneling amplitude in the Rabi frequency HQ = 



1/2 



2[(Ae(?j)/2) +v 2 (r])] ' . The last corrections scale as 
■q 2 at small values of 77, as it possible to see in Fig. [3] 
On the contrary, the energy difference increases linearly 
with respect to the asymmetry Ae(r?) ~ V((1r) — V(cll) + 
Hljr — fkuL = 2t)Vq + (3/8)r]huj. 

In the next example, we will consider a case in which 
the asymmetry affects substantially the tunneling dy- 
namics. 



B. Asymmetric parabolic potential 

The second example belongs to the class of asymmetric 
potentials with degenerate localized states el = £r. An 
example of such potentials is the following. For the neg- 
ative part (x < 0) we have V(x) = Vl(x) where Vl(x) 
reads 



V L {x) = V a 





\ 2 


(-- 


hi) -1 







, V 



(26) 



The energy of the localized states corresponds to El = 
—Vq + fajj/2. For the positive part x > 0, V(x) corre- 



sponds to V(x) — Vr(x) where Vn(x) reads 



V R (x) = V v 



±-1 



V v = 



(27) 



Similarly we have sr = —V n + hui^/2. This potential is 
continuous as Vl(0) = Vr(0) [39| . Generally we can have 
uj ri =/= uj and a v 7^ a, but the two levels are still degenerate 
under the condition 



1 



/ a v uj v \ 2 
V au; J 

A possible choice is 
w v =u(l-i]), a n /a = 



2Vo 



1 



0. 



(28) 



with the condition that i] < 2Vo/(hu). The result is 
shown in Fig. [4^,. 

Let us calculate now the tunneling amplitude 1/(77). As 
the degeneracy is not removed er = el, the prefactor A 
for the tunneling amplitude in Eq. ^ is one. In this case 
we have simply 



v(v) = vr = \/v(Q) v R {rj) 



(30) 



We have used = i/(0) as the left part of the potential 
is unmodified (see Fig. |4^i.). The tunneling amplitude 
associated to the left symmetric potential, Vl(x) = V(x) 
for x < 0, corresponds to the symmetric case with 77 = 0: 



V L = 1/(0) 




in 



(31) 



As the asymmetry does not change the shape of the po- 
tential, we can directly use the previous formula to obtain 
the tunneling amplitude associate to the right symmetric 
potential, Vr(x) = V{x) for x > 0, 




(32) 



The result of Eqs. f30 |l .(pT] ) . (p]l is shown in Fig. HJd, 
where we report again the ratio between ^(t7)/^(0) for 
different values of Vo/Hlj. We can observe that a small 
asymmetry of the potential 77 < 0.2 can rcnormalizc 
quantitatively the amplitude v of order 20%. As the 
Rabi frequency equals twice the tunneling amplitude 
HQ = 21/(77), m this case, the corrections due to the asym- 
metry appear clearly in the quantum dynamics of the 
two-level systems. 

The behavior of v shown in Fig. |4Jd can be explained 
by looking at the way we have chosen to deform the po- 
tential Eq. (|29|) . For the asymmetric potential, the bar- 
rier height at right is modified as [V(0) — V(<Xr)\/Vq = 
1 — rjHu / (2Vq) whereas the right attempt frequency is var- 
ied as ljr/uj =l — i]. The second correction dominates 
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Figure 4. (Color on line) a) The asymmetric parabolic po- 
tential Eqs. [|2S )) .l[27 |) .ip? ]) . b) The tunneling amplitude v(j]) 
as a function of the asymmetry parameter rj scaled with the 
symmetric value v(0) at different values of Vq/Tilj. 



owing to the condition hiu/(2Vo) < 1 for the two local- 
ized states. When 77 < 0, the right attempt frequency 
increases (hardening) leading to an enhancement of tun- 
neling. On the contrary, for 77 > the right attempt 
frequency decreases (softening) leading to a suppression 
of tunneling. A choice different from the one of Eq. (|29|) . 
for instance the one corresponding to increasing or de- 
creasing the height of the barrier in linear way as a func- 
tion of i], produces similar results but with a different 
dependence on r\. 



A time-dependent problem and quantum Zeno 
effect 



For the asymmetric parabolic potential with degen- 
erate local levels, we want to discuss a simple time- 
dependent problem. 

We assume to tune the asymmetry of the potential in 
time r\ = r](t) so that the tunnel amplitude gets a time 
dependence v = v(t). The time scale for the variations of 
r](t) (and therefore v(t)) are assumed slower enough to 
avoid excitations of the systems towards higher energy 
states of the potential E n with n > 2. Working under 
this assumption, the two-level description still holds. As 
a simple estimate, we consider the time scale for the vari- 
ations of v (t) smaller than the harmonic frequencies u>l 



and uir at the left and right well. Setting the energy 
level to £ = 0, we write the Schrodinger equation for the 
two-level system as 



dt { c R (t) 



-v{t) 
-v(t) 



CL(t) 

c R (t) 



(33) 



where Cf,(i) and c R (t) are the coefficients of the state 
at the time t. These equations can be easily solved by 
introducing the sum and the difference c± — Cl ± c R 
which satisfy the equation 



(34) 



Assuming that the system is prepared at the initial time 
t = in one localized state, let us say right, we have 

c L (t) = cos Q Jdt'v{t'fj , c R (t) = i sin Q Jdt'v(t'^J , 

(35) 

so that the probability to remain into the initial left state 
reads 



p L {t) = \c L {t)Y 



I - / dt'v{t') 



(36) 



An interesting case is the following evolution for v{t) 



V(t) =Vq + {Vl - Z/p) ^ Xn (t) , 



(37) 



where Xn{t) is the characteristic function equals to \n = 
1 in the time intervals [n(io+*i)+*o] < t < (n+l)(to+*i) 
and Xn = elsewhere, to is a time interval in which the 
tunneling amplitude is constant and equals to v = i>q 
whereas in a subsequent smaller interval t\ < to, the 
tunneling amplitude is suppressed to v\ <C vq [40| . The 
signal is repeated many times during a half Rabi period 
T/2 = n/Q = Trh/{2u Q ). 

The results for evolution of the probability Pr,(i) are 
shown in Fig. [5] both for the case v = v(£), Eq. (|37p . and 
for the case of constant amplitude v — vq. In Fig. (SJ 
we can observe that the probability Pz,[i, f(t)] is always 
higher then Pl[^, vq\. Reducing v to v\ during a short 
time interval t\ corresponds to deform the potential in 
time in a way to suppress the tunneling amplitude v\ <C 
vq. This corresponds to trap back the particle in the 
starting well at regular time intervals. As a consequence, 
a slow-down of the probability to escape from the initial 
well occurs, a result that is referred as quantum Zeno 
effect in the literature [29, K30I • 



V. CONCLUSION 

In summary, we have derived a useful and succinct 
expression for the tunneling amplitude v in asymmetric 
double- well potentials. Wc applied it to two examples: 
the quartic potential with a linear force and a kind of 
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t / (T/2) 

Figure 5. (Color on line) The probability Pt(t), full (red) line, 
Eqs. (|36l) . (|37|) for an asymmetric potential with degenerate 
local states. Parameters: vi/vq = 0.005, t\/to = 1/8 and 
to = T/16 where T = irh/(2v ). The dashed (black) line 
is the probability Pl calculated with constant amplitude vq. 
The dotted (blue) line represents the ratio u(t)/v(0). 

parabolic potential in which the asymmetry does not re- 
move the energy degeneracy of the two localized levels. 
From these simple examples, one can learn that there 
are no systematic effects of asymmetry on quantum tun- 
neling. The tunneling amplitude is enhanced or reduced 
depending on the shape of the potential and how the 
asymmetry is introduced. However we have illustrated as 
the formulas Eqs. ([I]),© allow one to obtain analytically 
and in a direct way the renormalization of the tunnel- 
ing amplitude in an asymmetric double-well potential in 
order to discuss its behavior as varying the asymmetry. 



ACKNOWLEDGMENTS 

The author thanks F.W.J. Hekking and R. Whitney for 
useful discussions and critical reading of the manuscript. 
This work was supported by ANR through contracts DY- 
COSMA and QUANTJO. G.R. acknowledges support 
from the European networks MIDAS, SOLID and GE- 
OMDISS. 



Appendix A: Comparison with exact numerical 
results 

In this appendix we show the comparison between the 
exact numerical result of the energy splitting (E\ — Eq) 
and the WKB semi-classical formula for the Rabi fire- 

1 /2 

quency Ml = 2[(e L (r}) ~ £ r(v)) + ^(v)} with v{rf) 
given by the Eqs. ©,© and the Eq. (fHj) . 

The first two eigenstates ipo (x) , (x) and their ener- 
gies Eq,Ei were computed numerically by discretizing 
the time-independent Schrodinger equation for the two 
potentials V(x) discussed in the paper. Using this ap- 
proach and imposing the boundary conditions ip n (x) = 
at the end points of a finite interval, the Schrodinger 



equation becomes a linear eigenvalue problem with a 
tridiagonal matrix assuming that the difference between 
the exact value of the second derivative d 2 t/j(x)/dx 2 and 
its discretized form [ip(xi+i) + tp(xi—i) — 2ip(xi)]/ Ax 2 is 
small [4l|. The first low-energy eigenstates are smooth 

1/2 

functions which extend over lengths cr s = [h/(mu) s )] 
Then one can choose a spacing Ax = min(p l , o r) / 'N 
with N sufficiently large (N ~ 10 3 ) in order to com- 
pute the eigenvalues with an acceptable error 8E\ , SE <C 
Ei — Eq. The end points were set to x m i n ~ —Mol on 
the left, and to x max = Man on the right (M ~ 5 — 10). 
A scaling analysis of Eq, Ei as functions of N, M was also 
carried out to test the convergence. 

In Fig. [S] we show two examples for the double- well 
parabolic potential with two degenerate levels £l = £r 
(Ml = 2v{rf}), i.e. Eqs. ([26]).(p7]l.(|29ll of Sec. ITVBl For 
the symmetric case r\ = in Fig. [5^,, the analytic semi- 
classical formula is in agreement with the numerical re- 
sult for Vo/(twj) > 0.4 (error < 1%) whereas for the 
asymmetric case rj = —0.5 in Fig. [6Jd the agreement is 
for V /(Mo) > 0.8. 



a) 




o [ — — t 

0.2 0.4 0.6 0.8 




Figure 6. (Color on line) Comparison between E\ — Eo nu- 
merically computed (the black dots with dashed line) and 2v 
(full red line) for the parabolic potential, a) Symmetric case 
r\ = 0. b) Asymmetric case r\ = —0.5. Inset: logarithmic scale 
for the y axis. 

In Fig. [7] we compare numerical and analytic solu- 
tions for the bias quartic potential Eq. ([23]) discussed 
in Sec. IIV Al In this case, the asymmetry due to the lin- 
ear bias r\ > removes the degeneracy of the localized 
states. As a consequence, at given rj > 0, the two-level 
approximation breaks down at large values of ratio Vq /fko 
as approaching the resonant condition el = £r + friOR 
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(see inset of Fig. Ub)- Thus the upper bound of valid- 
ity for our analytical approach is given by the condition 
V /(M < (cj R /u)Vo/(e L - e R ) ~ 1/(2??) to the leading 
order in -q and for Vo/(hiu) ^> 1. 
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Figure 7. (Color on line) Comparison between E\ — Eq nu- 
merically computed (the black dots with dashed line) and the 
semi-classical formula of the Rabi frequency (full red line) for 
the quartic potential, a) Case r\ — (Ml = 2v). Inset: loga- 
rithmic scale for the y axis, b) Case r\ = 0.2. Inset: a close-up 
for the range of large ratio Vo/hw. 
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